Ultrametric probe of the spin-glass state in a field 
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We study the ultrametric structure of phase space of one-dimensional Ising spin glasses with 
random power-law interaction in an external random field. Although in zero field the model 
in both the mean-field and non-mean-field universality classes shows an ultrametric signature 
[Phys. Rev. Lett. 102, 037207 (2009)], when a field is applied ultrametricity seems only present 
in the mean-field regime. The results for the non-mean field case in an external field agree with data 
for spin glasses studied within the Migdal-Kadanoff approximation. Our results therefore suggest 
that the spin-glass state might be fragile to external fields below the upper critical dimension. 
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I. INTRODUCTION 

Spin glasses^ 2 - are paradigmatic model systems that 
find wide applicability across disciplines. Although stud- 
ied intensely, our understanding of some of their funda- 
mental aspects is still in its infancy. In particular, the 
understanding of the nature of the spin-glass state re- 
mains controversial and active discussion has emerged 
recently.— ~— It is unclear if the mean-field replica sym- 
metry breaking (RSB) picture^ of Parisi describes the 
non-mean-field behavior of spin-glasses in an externally- 
applied field best. While the droplet theoryi^r— states 
that there is no spin-glass state in a field for short-range 
systems, the mean-field RSB picture 2 ^ 14 i 19 i 20 states that 
for low enough temperatures T and fields H (i.e., below 
the de Almeida-Thouless line)2i a stable spin-glass state 
emerges. The question lies at the core of theoretical de- 
scriptions and is of immediate importance to applications 
in research fields ranging from, for example, sociology to 
economics where terms linear in the spin variable can 
emerge. 

One way to settle the applicability of the RSB picture 
to short-range spin glasses in a field while avoiding tech- 
nical difficulties when measuring observables in a field ^ 
is by testing^ 2 - if the phase space is ultrametric (UM). 
Unfortunately, the existence of an UM phase structure 
for short-range spin glasses on hypercubic lattices re- 
mains elusive^ mainly because only small systems can 
be studied numerically. Recent results in zero field 2 ^ sug- 
gest that short-range systems are not UM, whereas other 
opinions exist i 2 ^— 

More recently^ 2 - results on one-dimensional (ID) Ising 
models with power-law interactions showed that short- 
range spin glasses might be UM after all. Therefore, a 
natural probe for the spin-glass state in a field is to study 
the UM response of ID Ising models with power-law in- 
teractions when an external field is applied. The model 
has the advantage in that by tuning the exponent of the 
power law, the universality class can be tuned between 
a mean-field and a non-mean-field regime. In addition, 



large linear system sizes can be simulated, which allows 
for a better finite-size scaling analysis than for hypercu- 
bic lattices^ 

Our results show that for this model in a field the phase 
space has an UM structure in the mean-field regime. 
However, in the non-mean-field regime, when an exter- 
nal field is applied, the UM structure seems to be much 
weaker for the studied system sizes, suggesting that the 
spin-glass state for short-range systems is fragile with 
respect to externally-applied fields. These results are 
compared to studies of spin glasses within the Migdal- 
Kadanoff (MK) approximation. 




FIG. 1: Dendrogram obtained by clustering 100 con- 
figurations (see text) for a sample system with a — 0.0 
(Sherrington-Kirkpatrick model) and L — 512 at T = 0.36, 
together with the matrix d a p (grayscale, distance is black) . 
The order of the states is given by the leaves of the dendro- 
gram (figure rotated clockwise by 90°). 



II. MODEL 

The ID Ising chain with long-range power-law 
interaction o 17 i 27 ~ — is described by the Hamiltonian 



H 



^ t JijSiSj ^ ] hiSi ; Jij — c(u)- 



(1) 



2 



where Si G {±1} arc Ising spins and the sum ranges 
over all spins in the system. The L spins are placed 
on a ring to ensure periodic boundary conditions and 
Tij = (L/tt) sin(7r|i — j\/L) is the geometric distance be- 
tween the spins, are Gaussian random couplings. The 
constant c(cr) is chosen^ such that for the mean-field 
transition temperature T c MF (cr < 0.5, L,H = 0) = 1. In 
Eq. ([T]), the spins couple to site-dependent random fields 
hi chosen from a Gaussian distribution with zero mean 
and standard deviation [hf}H 2 = H. 

The model has a rich phase diagram when the expo- 
nent a is changed^ both the universality class and the 
range of the interactions can be continuously tuned. In 
particular, a = gives the Sherrington-Kirkpatrick (SK) 
model, 3 - whose solution is the mean-field theory for spin 
glasses and where a spin-glass state in a field is expected 
(i.e., an UM signature for low enough H and tempera- 
tures T). More importantly^ for 1/2 < a < 2/3 the 
critical behavior is mean- field-like, while for 2/3 < a < 1 
it is non-mean-ficld-likc. 

Here we study in a field H = 0.10 the SK model 
[a = 0] to test our analysis protocol, as well as the 
ID chain for a = 0.60 (also mean- field-like) , as well as 
a = 0.75 (T c r~i 0.69, roughly corresponding to four space 
dimensions) outside the mean-field regime. We choose 
two values of a ^ to be able to discern any trends 
when the effective dimensionality 3 ^ is reduced. In general 
d e s = (2— r/)/(2<7— 1), where r\ is the critical exponent for 
the short-range model at space dimension d = d e g. Note 
that r\ is zero in the mean-field regime and, for example, 
-0.275(25) for d = 4M 



III. NUMERICAL METHOD AND 
EQUILIBRATION 

We generate spin-glass configurations by first equili- 
brating the system at low temperatures and an exter- 
nal random field of standard deviation if = 0.1 using 
the parallel tempering Monte Carlo methodi 33 i 34 Once 
the system is equilibrated we record states ensuring that 
these are well separated in the Markov process and thus 
not correlated. In practice, if we equilibrate the system 
for -r cq Monte Carlo sweeps, we generate for each dis- 
order realization 10 3 states separated by r cq /10 Monte 
Carlo sweeps. We test equilibration using the method 
presented in Ref. [ill We consider systems sizes up 
to L = 512, which is the same maximum size as in 
the zero-field case studied previously^ but numerically 
much harder than in the zero-field case because Monte 
Carlo methods equilibrate considerably slower in a field. 
For the parallel tempering simulations T m j n = 0.36 and 
Tmax = 1-40 (16 temperatures). For all values of a stud- 
ied, and all system sizes L, we generate 4000 disorder 
realizations. For L = 32, the equilibration time is 2 x 10 4 
Monte Carlo sweeps (MCS), for 64, 1.5 x 10 s MCS, for 
128, 5 x 10 5 , and for 256 and 512, 10 6 MCS. 

The presented data are for T = 0.36. In Ref. HH we 



fixed T ss 0.4T C for all values of a studied to ensure that 
we are deep in the spin-glass phase. However, it is unclear 
if one-dimensional spin glasses with power-law interac- 
tions have a spin-glass state in a field for a > 2/3 i 11 ' 13 i 36 
Using the T c estimates of Leuzzi et al^ at zero and finite 
field (H — 0.1) for the diluted version of the model we 
estimate that if a spin- glass state exists for H = 0.1 it 
should suppress the zero-field T c by approximately 20%. 
For a = 0.75 it is known that T C (H = 0) « 0.69(l).ii 
Therefore T = 0.36 corresponds roughly to a 40% reduc- 
tion of the critical temperature (i.e., deep in the putative 
spin-glass phase). 

We also study spin glasses within the standard 
MK approximation^ (i.e., spin glasses on hierarchical 
lattices) Due to the simple lattice structure, the phase 
space is also expected to be simple. In fact, as shown rig- 
orously in Ref. [39|, spin glasses on MK lattices are replica 
symmetric. We used a variation of the standard MK re- 
cursion where, starting from one bond, iteratively each 
bond is replaced by 2 d bonds and 2 d ~ 1 spins (d = 3). For 
details, see, for example, Refs. and EH. 



IV. ULTRAMETRICITY 

Ultrametricity appears in different fields of research 
ranging from linguistics to the taxonomy of animal 
species and is a key component of Parisi's mean-field so- 
lution of the SK modeli^ ' 14 i 42 Therefore, if a spin glass 
has no UM phase-space structure there is a strong indi- 
cation that Parisi's mean-field picture might not work for 
this system. 

In an UM spaced 3 - the triangle inequality d ai < d a p + 
dpy is replaced by a stronger condition where d ai < 
max{d a/ 3, d^} (i.e., the two longer distances must be 
equal and the states lie on an isosceles triangle). Here, 
d a p represents the distance between two points a and ft 
in phase space. 

We use the approach developed in Ref. [TU which 
is closely related to the one used by Hcd et al. in 
Ref. HH. For each disorder realization we produce M = 
10 3 equilibrium configurations. These are sorted using 
the average-linkage agglomerative clustering algorithm! 4 ^ 
The clustering procedure starts with M clusters contain- 
ing each exactly one configuration. Distances are mea- 
sured in terms of the Hamming distance d a p = (1 — \q a fs I ) , 
where q a p = N^ 1 J^i S?S? is the spin overlap between 
configurations {S a } and {S"}. Iteratively the two closest 
clusters C a and C& are merged into one cluster d, reduc- 
ing the number of clusters by one. The distances of the 
new cluster Cd to the other remaining clusters have to be 
calculated: The distance between two clusters is the aver- 
age distance between all pairs of members of the clusters. 
The iterative procedure stops when only one cluster re- 
mains, the results are then typically structured in a tree- 
like structure called a dendrogram (see Fig. [lj . To probe 
for a putative UM space structure, we randomly select 
three configurations from the hierarchical cluster struc- 
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FIG. 2: (Color online) Distribution P(K) for different system 
sizes (all panels have the same horizontal and vertical scale) 
and an external random field H = 0.1. (a) Data for the SK 
model. The distribution diverges very slightly for K — > 
and L — > oo thus signaling an UM phase structure, (b) Data 
for a = 0.60 (mean-field universality class). There is still a 
weak hint of a divergence for K — > 0. (c) Data for a = 0.75 
(non-mean-field universality class) . There is no clear sign of a 
divergence in P(K) for K — > 0. Note that when H — data 
for a = 0.75 show a clear signature for UM behavior.— Error 
bars are smaller then the symbol size. 



ture (see Ref. l22f ). resulting in three mutual distances. 
Next, we sort these Hamming distances d max > d mo d > 
d m i n and compute K — (d max — d mc d) / g(d) , where g(d) is 
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FIG. 3: (Color online) Variance Va,r(K) of P(K) as a function 
of system size L for different values of a. The data can be 
fit to a power law (dashed lines). In the mean-field regime 
(SK and a — 0.6) a fit to a constant is unlikely (see text). 
The power-law decay of the variance as a function of system 
size suggests a divergence in P(K) for K — > 0. For a — 0.75 
the data are compatible with a constant (solid line) or a very 
weak power-law behavior. 



the width of the distance distribution. If the phase space 
is UM, then we expect d max = dmcd for L — > oo. Thus 
P(K) — > S(K = 0) for L — >• oo and the for the variance 
of the distribution Vax(K) — > for L — > oo. 
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FIG. 4: (Color online) Variance Var(Tf) of P(K) as a function 
of system size L for spin glasses on MK lattices. The data are 
compatible with a constant behavior, showing that there is 
no UM phase-space structure for spin glasses within the MK 
approximation. The solid line is a guide to the eye. 



V. RESULTS 

Figure EJa) shows the distribution P(K) for the SK 
model (cr = 0), T = 0.36, and H = 0.10. There is a slight 
hint for a divergence for K — > 0. Similar results are found 
for the mean-field regime with a = 0.60 [Figure (2jb)]. 
The UM signature in a field is considerably weaker than 
when no field is applicdfi^ While for the SK model there 
is still a faint sign of a divergence, for larger values of a it 
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is hard to see if the distributions diverge for K — > and 
L — > oo. Figure [He) shows data for cr = 0.75, T = 0.36, 
and -ff = 0.10 where no clear sign of a divergence is 
present, suggesting that phase space might not be UM 
outside the mean-field regime. 

Hence, drawing conclusions from the P(K) data is not 
sufficient. A better probe is given by the variance Var(A") 
of P{K) as a function of system size L (Fig. The 
variance of the distribution for the SK model clearly de- 
cays with a power law Var(A) ~ [6 = 0.49(4), 
7 = 0.13(2), Q-factor - 0.28] 4^ If we restrict the fit to 
L > 128 we obtain b = 0.58(7) and 7 = 0.16(2) with a Q- 
factor ~ 0.487. A fit to a constant gives Q = if the fit is 
performed for all data or restricted to L > 128. A fit to 
a constant+power-law behavior Var(A) ~ a + b/L" / gives 
a constant a compatible with zero and a clear power-law 
decay. Therefore, and as expected, the SK model shows 
an ultrametric phase space structure for small externally 
applied magnetic fields. 

Similar results are obtained for a = 0.60 where a fit 
to a power law is very likely with b — 0.395(6), 7 = 
0.074(3), and Q = 0.989 [restricted to L > 128 we obtain 
b = 0.374(1), 7 = 0.064(1), and Q = 0.983]. However, 
a fit to a constant gives Q < 10 -5 (0.124 restricted to 
L > 128). We also attempted a fit to a constant+power- 
law behavior [i.e., Var(A') ~ a + b/ L 1 ]. We obtain a = 
0.18(2) > with Q = 0.989. This suggests that we might 
be at a marginal regime (i.e., close to the upper critical 
dimension) . 

For a = 0.75 a fit to a very weak power law with b = 
0.30(1) and 7 = 0.014(6) is found with Q = 0.897. Thus, 
the exponent 7 is extremely small, only within about two 
standard deviations from zero. Correspondingly, a fit to 
a constant is equally probable with Q = 0.811. Similar 
results are obtained for L > 128 where b = 0.33(2) and 
7 = 0.028(9) with Q = 0.811, and Q = 0.766 for a fit to a 
constant. A fit to a constant+powcr-law behavior gives 
a power-law exponent consistent with zero within error 
bars. 

Summarizing, either ultramctricity in the non-mean- 
field regime is completely lost in a field or greatly weak- 
ened, suggesting a marginal signal for a = 0.60. Larger 
systems would be needed to fully discern the behavior, 
however they are out of reach with current technology. 
Note that for diluted systems larger system sizes are pos- 
sible, but the finite-size effects arc stronger, resulting in 
no overall benefit . 



Within the MK approximation the distributions P(K) 
also show no divergence for K — > 0. Figure |4] shows the 
variance of the distributions as a function of the system 
size for very large lattices. There is no discernible de- 
crease with an increasing number of spins (i.e., no UM 
structure of phase space). In fact, a fit to a power-law 
behavior results in a slope compatible with zero (i.e., a 
constant behavior). This is to be expected because the 
model is defined on a hierarchical lattice. However, a 
direct comparison to the results for a = 0.75 strength- 
ens the evidence of a potential non-UM structure for the 
latter case, in agreement with recent results^ 



VI. SUMMARY AND CONCLUSION 

We have studied numerically the low-temperature 
configuration landscape of long-range spin-glasses with 
power-law interactions. By tuning the exponent a that 
governs the decay of the power-law interactions and 
therefore their range we can tune the system out of the 
mean-filed universality class. Using a hierarchical clus- 
tering method and analyzing the resulting distance ma- 
trices we show that when a field is applied the system 
is only clearly UM in the mean-field regime, unlike in 
the zero-field case where an UM signal was found for 
values of a that correspond to space dimensions above 
and below the upper critical dimension. Therefore, our 
results suggest that the spin-glass state is fragile to an 
externally-applied field below the upper critical dimen- 
sion. Larger systems would be needed to determine if 
the UM signature for a = 0.75 (corresponding approx- 
imately to four space dimensions) persists in a field or 
not. 
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